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Abstract. We present a spectroscopic analysis of a sample of 21 galaxies with z = 0.2-1.5 drawn from a 25 
square arcmin ultra-deep ISOCAM survey at \ e ff = 15 fJ.m centered in the WFPC-2 Hubble Deep Field South. 
• Near-infrared spectra are reported for 18 ISO sources, carried out with ISAAC on the VLT, aimed at detecting 

the redshifted H Q +[Nll]. Additional optical data come from the ESO VLT/FORS2 and NTT/EMMI, primarily 
q . targeting [On], [Om] and H^ for further physical insight. Although not numerous in terms of areal density in 

the sky, this population of very luminous IR sources has been recently found to be responsible for a substantial 
fraction of the extragalactic background light energy density. H a line emission is detected in virtually all the 
observed objects down to a flux limit of 7xl0 -17 erg cm -2 s _1 (corresponding to LH a > 10 41 erg s _1 at z = 0.6 
for H = 65, £Ia = 0.7 and Q m = 0.3). Our analysis (including emission line, morphology, and SED properties) 
shows clear evidence for AGN activity in only two of these sources: one type-I (with broadened H Q at z=1.57) 
and one type-II quasars (with inverted [Nll]/H a ratio at z=1.39), while we suspect the presence of an AGN in 
two further sources (an Ultra-Luminous IR Galaxy, ULIRG, at z=1.27 and a luminous galaxy at z=0.69). The H a 
luminosities indicate star formation rates (SFR) in the remaining sources between 0.5 and 20 Mg/j/r, assuming a 
Salpeter IMF between 0.1 and 100 Mq and without extinction corrections. We find good correlations between the 
mid-IR, the radio and H a luminosities, confirming the mid-IR light as a good tracer of star formation (while the 
SFR based on H a flux show some large scatter and offset, which are still to be understood). We have estimated 
the baryonic masses in stars with a newly-developed tool fitting the overall optical-IR continuum, and found that 
the host galaxies of ISO sources are massive members of groups with typically high rates of SF (SFR ~ 10 to 300 
Mq/uv). We have finally compared this ongoing SF activity with the already formed stellar masses to estimate 
the timescales tsF for the stellar build-up, which turn-out to be widely spread in these objects between 0.1 Gyrs 
to more than 10 Gyr. The faint ISOCAM galaxies appear to form a composite population, including moderately 
active but very massive spiral-like galaxies, and very luminous ongoing starbursts, in a continuous sequence. From 
the observed tsF and assuming typical starburst timescales, we infer that, with few exceptions, only a fraction of 
the galactic stars can be formed in any single starburst event, while several of such episodes during a protracted 
SF history are required for the whole galactic build-up. 
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1. INTRODUCTION 

New populations of faint high-redshift sources have been 
recently discovered by deep mid- and far-IR surveys with 
the Infrared Space Observatory (e.g. Genzel & Cesarsky 
2000, Aussel et al. 1999, Altieri et al. 1999, Elbaz et al. 
1999, Rowan-Robinson ct al. 1997) and by large millimet- 
ric telescopes at longer wavelengths (e.g. Smail et al. 1997, 
Hughes et al. 1998). 

These sources display various distinct features com- 
pared with other optically selected galaxy populations. 
They are very luminous on average (Ltoi > 10 11 Lq, Elbaz 
et al. 2002), with the bulk of their emission coming out in 
the far-IR, in a similar way as the IRAS-selected galaxies 
include the most luminous systems in the local universe. 
On the contrary, their areal density (a few sources/square 
arcmin at the faintest limits detectable by ISO) is much 
lower than found for faint blue galaxies in the optical (e.g. 
Ellis 1997). 

Another remarkable property of the faint IR-selected 
sources is to display extremely high rates of evolution with 
redshift, exceeding those measured for galaxies at other 
wavelengths and comparable or larger than the evolution 
rates observed for quasars (Hughes et al. 1998, Barger 
et al. 1998, Elbaz et al 1999, 2002, Blain et al. 1999, 
Franceschini et al. 2001). This fast evolution of the IR 
sources implies that dust-obscuration, if moderately im- 
portant in local galaxies where less than 50% on average 
of the optical-UV emission is absorbed, has strongly af- 
fected instead the past active phases of galaxy evolution. 

Franceschini et al. (2001) and Elbaz et al. (2002) have 
matched the statistical and IR-spectral properties of the 
faint ISO sources detected at A e // = 15 /im with the 
spectral intensity of the recently discovered Cosmic IR 
Background (CIRB, see Hauser et al. 1998, Puget et al. 
1996). The CIRB appears to contain a large fraction (up 
to ~ 70%, though this number is made uncertain by 
the optical-UV background intensity) of the total extra- 
galactic background energy density from radio to X-rays, 
Cosmic Microwave Background excluded. These analyses 
have found that, due to their high luminosities and mod- 
erate redshifts (z ~ 0.5 to 1.3), the faint 15 fim sources 
include the main contributors to the CIRB. This is a ro- 
bust conclusion, based on the observed shape of the 15 /im 
counts, see Elbaz et al. 1999, and only assuming for these 
sources a typical IR galaxy SED. Then a large fraction 
of stars in present-day galaxies, or alternatively the bulk 
of degenerate baryons contained in nuclear supermassive 
BH's, have formed during IR-luminous dust-cxtinguished 
evolutionary phases. 

The deep diffraction- limited observations at 15 yum 
with ISO provide quite an effective way of probing this 

* Based on observations with ISO, an ESA project with in- 
struments funded by ESA member states (especially the PI 
countries: France, Germany, the Netherlands, and the United 
Kingdom) with the participation of ISAS and NASA. 
** Based on observations collected at the European Southern 
Observatory, Chile, ESO No 63.O-0022, 65-000 67-000 



dust-obscured high-redshift activity. The numerous source 
samples detected in this way offer an important advantage 
over longer wavelength observations to allow easy opti- 
cal identification, thanks to the relatively small error-box 
(4.6" PSF, Okumura 1998) and the moderate redshifts 
and faintness of the optical counterparts. 

We report here on optical and near-IR spectroscopic 
follow-up of a representative subset of the faint ISO pop- 
ulation selected from a region centered in the HDF South. 
Apart from measuring the redshift, motivation for our ob- 
servations was to clarify the nature of these objects: the 
main open issues were to assess the presence of energet- 
ically dominant AGNs as power sources, and to estimate 
the main physical parameters of the starburst and normal 
galaxy population, like the Star- Formation Rate (SFR), 
the extinction, and the stellar mass. Section 2 describes 
the IR-sclected sample and the observations. Section 3 
analyses the source properties based on the emission lines, 
while Section 4 those of the SEDs and continuum emis- 
sion. Section 5 compares optical, near-IR, far-IR and ra- 
dio indicators of SF and discusses the level of activity in 
the IR-selected galaxies. Section 6 summarizes our conclu- 
sions. 

All quantities are computed assuming a universal ge- 
ometry with H = 65 km s _1 Mpc -1 , £l m = 0.3, £!a = 0.7. 
We indicate with the symbol S15 the flux density in Jy at 
15 \xm (and similarly for other wavelengths). 

2. THE SAMPLE AND THE OBSERVATIONS 

2.1. Sample selection 

The Hubble Deep Field South was observed between 
October 17 and November 29, 1997, with the array cam- 
era ISOCAM onboard the Infrared Space Observatory, as 
part of the ELAIS collaboration (Oliver et al. 2000). The 
observations were carried out with two broad-band filters, 
LW2 (5-8.5 /um, X eff = 6.75 /xm) and LW3 (12-18 /urn, 
X e ff = 15 ^m). This followed a previous similar observ- 
ing campaign on the HDF-North (Oliver et al. 1997), but 
adopted an improved observing strategy, particularly with 
LW2. The deep ISOCAM images with the two filters, ob- 
tained as repeated raster scans to improve the flat-field ac- 
curacy and time-redundancy, covered the same area cen- 
tered on the HST WFPC-2 field. All details on the ob- 
servations can be found in Oliver et al. (2002). We will 
consider in the following only the LW3 sample selected at 
15 //m (for these sources we will also make use of 6.75 /xm 
fluxes or upper limits from the LW2 observation). 

Oliver ct al. (2002) and Aussel et al. (2003) analyzed 
the ISOCAM data with two independent methods, ac- 
counting in detail for the time-varying signals under the 
effect of cosmic ray impacts. In particular, the method 
adopted by Aussel ct al. makes use of a wavelet analysis 
of the combined spatial-temporal observable space (the 
PRETI method, see Stark et al. 1999, Aussel et al. 1999). 
The PRETI reduction has detected with LW3 63 sources 
brighter than Si5 Mm = 90 fiJy (59 above 100 fiJy) over an 
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Fig. 1. Image of the HDF-S region with indicated the LW3 15 /xm sources (source labels refer to the catalogue by 
Aussel et al. 2003, see Tables Q] and following). Of the 86 objects indicated in this map, 63 belong to the complete 
sample with Sis > 90/iJy, while the other are fainter than this limit. The map is a collage bv lHoolj l|l999|) including 
WFPC-2 F814 images in the Flanking Fields. The map scale is 4.7 arcminutes on a side. Contours represent the LW3 
image depth, increasing towards the center. 



area of 25 square arcminutes (for comparison, 24 of these 
sources appear in the shallower list by Oliver et al. within 
the inner 19.6 square arcmins). Detailed simulations have 
shown that above Si5 M , n = 100 /xJy the PRETI sample is 
complete and free of spurious sources. 

2.2. Source identification, target selection, 
photometric redshifts 

We have compared the ISOCAM source lists with those 
from the Deep ESO Imaging Survey (EIS Deep), includ- 
ing optical imaging in UBVRI with ESO NTT/SUSI-2 
down to limiting magnitudes of Vab ~ 27, Bab ~ 26.5, 



V ab ~ 26, Kab ~ 26, Iab ~ 25 and near-infrared JHK 
observations down to Jab ~ 25, Bab ~ 24, Kab ~ 24 
performed with NTT/SOFI (Da Costa et al. 1998). In 
ad dition we have used the optical catalogues in uBVRI 
bv iTeolitz et a l. (1998), which has 5er limiting magni- 
tudes of u~24.5, B~26.1, V~25.5, R-25.4 and 1-24. 
Unfortunately, the EIS Deep images cover only a fraction 
(~70%) of the ISOCAM survey area. 
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Table 1. Available photometric data for the 15 /im sources in the HDF-S observed spectroscopically in the optical and near- infrared. Data are from the EIS 
Deep survey, but in highlighted cases for which the Teplitz et al. (1998) photometry is used. All magnitudes are in the AB system; ISOCAM LW2 and LW3 
fluxes and uncertainties are reported in the last four columns. In the case only an upper limit is available, it is specified. 



Franceschini, Berta, Rigopoulou, et al.: IR Spectroscopy of Faint ISO Galaxies 



5 



The spatial resolution of the ISOCAM images corre- 
sponds to a PSF of 4.6" (Okumura 1998). Within the 
ISOCAM errorbox, it turns out that there is almost in- 
variably a galaxy relatively bright in the red wavebands 
(I,J,H,K). The procedures for optical identification are de- 
tailed in Mann et al. (2002) (see also Aussel et al. 1999). 

Of the 63 sources in the PRETI LW3 sample, 24 are 
outside the JHK EIS coverage (10 of these are also outside 
the EIS optical imaging). 

Of the other 39 sources, four are galactic stars, while 
for the remaining 35 objects the optical/near-IR cover- 
age of the galaxy SED is detailed enough to allow a pre- 
cise estimate of the redshift through fits with spectro- 
photometric models. 

We selected the targets for the VLT/ISAAC spectro- 
scopic follow-up from the HDF South LW3 source list 
based on the following criteria: a) H a should be in the 
wavelength range covered by the ISAAC gratings, b) a 
secure counterpart should exist in the I or the K band 
images. 

We did not apply any selection based on redshifts or 
colours, except to ensure that H Q was within the ISAAC 
spectral range. With these constraints, our reference sam- 
ple reduces to 25 galaxies with 15 /im flux densities rang- 
ing between 95 and 800 /j,Jy. It is thus a representative 
sample of the strongly evolving ISOCAM population near 
the peak of the differential source counts (Elbaz et al. 
1999). 

From these 25 sources we randomly selected 18 for the 
ISAAC follow up. For part of these and for 3 additional 
objects we have optical spectroscopic data. 

The source list, coordinates, and the available opti- 
cal, near-IR and mid-IR photometric data are reported in 
Table □ Stamp images from the F814 WFPC-2 maps for 
each sources are reported in Appendix together with 
plots of the optical-IR SEDs and notes on the individual 
sources. 

To set the ISAAC grating (Z, SZ, J, H) for H a de- 
tections in our target sources, we used redshifts from op- 
tical spectra available for some of the objects at z<0.7 
(Dcnnefeld, private communication, see also Rigopoulou 
et al. 2000). 

For all other sources we used photometric redshift es- 
timates based on fits of the observed SEDs with synthetic 
galaxy spectra. In particular, we have used PEGASE (Fioc 
and Rocca-Volmerange 1997) and GRASIL (Silva et al. 
1998) to construct grids of spectra as a function of galactic 
age, considering two evolutionary sequences: one describ- 
ing a spiral- like evolutionary model with exponentially de- 
creasing star-formation and long (5 Gyr) timescale, the 
other reproducing a more typical evolution for ellipti- 
cals with short (1 Gyr) exponential timescale and star- 
formation truncated by a galactic wind. 

Since the main spectral feature of the SED suited for 
the redshift estimates is the Balmer 4000 A discontinuity, 
the effect of dust extinction on the redshift estimate is 
modest. The fit is automatically found from x 2 interpola- 
tion on a 2D grid of redshift versus galactic age. 



HDFS: LW3 sources 



Fig. 2. The redshift distribution of faint ISOCAM LW3 
sources in HDF South with optical identifications. The 
objects range between z=0.2 and 1.6, mostly due to K- 
correction effects. Redshifts are spectroscopic when avail- 
able, otherwise photometric. 
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*: high-res data (R a ~ 5000) also exist. 
Table 2. Summary of ISAAC/VLT observations. Slits 
have been rotated (anti-clockwise with respect to the 
North direction by the angle in column PA) in order to 
include the most interesting features of the sources. 



After spectroscopic confirmation, our photometric red- 
shift determinations turned out to be accurate to Az~ 
±0.1. The fits based on the GRASIL and PEGASE codes 
provided essentially the same results in terms of redshifts. 
The distribution of redshifts for the ISOCAM LW3 sam- 
ple in the HDF South is shown in Fig. |2 The objects 
range between z=0.2 and 1.6, an interval mostly imposed 
by K-correction effects (e.g. Franceschini et al. 2001): at 
z ~ 0.6 — 0.8 the 7.7 /im PAH emission feature usually 
very intense in starburst galaxies falls within the LW3 fil- 
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tor. The z distribution shows a strong peak at z ~ 0.6 
probably due to a cluster or a large galaxy concentration 
in the HDF South. This overdensity is also apparent in 
the analysis of photometric redshifts by Rudnick et al. 
(2001). This z distribution is remarkably different from 
that of ISOCAM sources at similar depths in the HDF 
North (Aussel et al. 1999), particularly in lacking the peak 
at z~ 0.8 — 1. The effect of cosmic variance is very promi- 
nent between the two fields. 

2.3. IR spectroscopic observations and data reduction 

We collected the IR spectra during three runs (September 
1999, August 2000 and 2001) using the infrared spectro- 
graph ISAAC (Moorwood et al. 1998) on the ESO ANTU 
telescope (formerly UT1), on Paranal, Chile. The observ- 
ing logs are sumarized in Table El 

Observations have been performed with the Low 
Resolution grating, providing a spectral resolution 
R s -600 for a slit width of 1" (the length is fixed to 2'). 
Four of our sample objects (S55, S27, and an interacting 
pair associated with source S25) were also observed dur- 
ing the 2000 and 2001 runs with the Medium- Resolution 
gratings (R s -5000) (Rigopoulou et al. 2002). 

To maximize the observing efficiency, whenever possi- 
ble the slit position included two target galaxies at any 
given orientation. Most of the targets were first acquired 
directly from a 1-2 min exposure in the H-band. In the 
case of the very faint objects (H > 20.0 mag), we blind- 
offset from a brighter star in the HDF-S field. Observations 
were made by nodding the telescope ± 20" along the slit 
to facilitate sky subtraction (always avoiding overlap of 
the two objects in the slit). Individual exposures range 
from 120 (in H and J bands) to 300 (SZ band) seconds. 
During the 1999 and 2000 runs, sky conditions were excel- 
lent throughout the acquisition of the spectra, with see- 
ing values typically in the range 0.4"-0.8 and dipping 
down to 0.25". For each filter, observations of spectro- 
scopic standard stars were made in order to flux calibrate 
the galaxy spectra. 

The data were reduced using applications in the 
ECLIPSE (Devillard 1998) and IRAF 1 packages. Accurate 
sky subtraction is critical to the detection of faint lines. 
Sky was removed by subtracting the pairs of offset frames. 
In some cases this left a residual signal (due to tempo- 
ral sky changes) which was then removed by performing 
a polynomial interpolation along the slit. OH sky emis- 
sion lines were also carefully removed from the spectra. 
Spectrum extraction for each galaxy was performed using 
the APEXTRACT package. Standard wavelength calibra- 
tion was applied. 

The ISAAC spectra have bee n flux-ca l ibrate d us- 
ing standard infrared stars from IPicklesI l|l998l) and 



1 The package IRAF is distributed by the National Optical 
Astronomy Observatory which is operatedby the Association of 
Universities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Fundation. 



Ivan der Bliek et"afl l|l996|) . The near- infrared spectra 
from the 1999 run have appeared in Rigopoulou et al. 
(2000). Here we present in Figure the ISAAC spectra 
from the 2000 and 2001 runs, while the measured fluxes 
are summarized in Table [3] 

2.4. Complementary optical spectroscopy 

To add further constraints on the nature of the 15 
/jm source population, particularly to estimate dust- 
extinction from Balmer line ratios, two different sets 
of ground-based optical spectroscopic observations have 
be en analyzed. 

iTresse et all dl999|) observed the region around QSO 
J2233-606, which is in the STIS HDF-South, with the 
ESO Multi Mode Instrument (EMMI, D'Odorico 1990) at 
NTT, during the nights between September 23 th and 25 th 
and October 17 th to 19*' 1 1998. The data were obtained 
in Multi Object Spectroscopy (MOS) mode, with different 
pointings and masks, for a total of 6 sets and roughly 200 
slit positions. ITresse et alJ l)l999|) chose slits with 1.02" or 
1.34" width and spectral resolutions of 10.6 and 13.9 A in 
the I band (centered at 7985A). 

We also obtained observations with FORS2 
(Appenzeller et al. 2000) at VLT/UT2 in MOS con- 
figuration, during the FORS2 commissioning phase, 
between December 22, 1999, and January 5, 2000. The 
300/ grism with a dispersion of 2.5A/pixel between 6000 
and 11000 A and a l"-wide slit were used. A few of the 
ISOCAM targets in our list have been included by chance 
in this spectroscopic survey. 

Both EMMI and FORS2 data were reduced with the 
standard IRAF's tasks. Pre-reduction (bias subtraction 
and flat -fielding) , extraction and wavelength calibration 
were performed with the usual proce dure. EMMI data 
were flux- calibrated using stars f rom IStone fc Baldwinl 
l|l983h and lBaldwin fc Stonel l|l984h . observed within each 
mask. Unfortunately, for FORS2 observations no data on 
standard stars were available. However, two stars have 
been included in the MOS HDF-South frames: after iden- 
tification in the EIS catalogue, the spectral data were 
calibrated by imposing to the observed stellar spectra 
to reproduce the I-band fluxes, after convolution with 
the instrument response function reported in the ESO 
exposure-time calculator. 

Optical spectra with EMMI are reported in Figure 0] 

3. SOURCE PROPERTIES BASED ON THE 
EMISSION LINES 

The ISAAC near-IR spectroscopy has been targeted to 
detect the redshifted H a line, while optical observations 
with EMMI and FORS2 have allowed the study of the 
rest-frame emission lines up to ~ 5000A. 

Table summarizes the properties of reliable line de- 
tections with the three instruments. The columns report, 
in the order, the source identification, the instrument, in- 
tegration times, the measured fluxes of [On]A3727, Hg, 
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: data from Ris'ODOulgu^^il^ (2000) 



Table 3. Summary of results of spectroscopic observations. We report the source name, the instrument, the exposure 
time, spectroscopic redshift, the measured fluxes (in units of 10~ 17 erg cm~ 2 s _1 ) and equivalent widths (EW, in 
A) of [Oil1A3727 R0 ( A 4861), [Om]AA4959, 5007, and Ha (A 6563). Ha fluxes and EW of Ha+[Nn] are from 
IRigopoulou et all l|2000|) . These measured line fluxes are not corrected for aperture. 



[Om]AA4949, 500 7, H„+[Nn1. The T able i ncludes revised 
fluxes from the iRigopoulou et al"! l)2000j) observations. 
However, only the extremely good seeing conditions of 
the 1999 observing run have allowed the separation of the 
[Nil] from H a for some of the sources. The correspond- 
ing separate fluxes for t he resolved lines are reported in 
IRigopoulou et all (2000). If an object was detected in dif- 
ferent runs, both results are shown. Uncertainties on the 
measured fluxes are of the order of 10%. 



3.1. Aperture corrections 

We have applied aperture correction factors to account for 
the light missed due to our finite 1" size slit. We have cal- 
culated different correction factors for unresolved (source 
sizes < 1") and resolved sources (size > 1"). 

For the unresolved sources we applied an average cor- 
rection factor of 1.15 which was calculated as follows. We 
have smoothed the WFPC2 F814 image of the galaxies 



to a resolution of 0.6" (to simulate the seeing during the 
observations) and measured the ratio between the total 
source counts and the counts through an artificial 1" slit. 

For the resolved sources (S27, S28, S53 ands S55) aper- 
ture corrections have been calculated with two different 
methods. 

(1) We used the WFPC2 F814 image (smoothed as 
before to a resolution of 0.6" to account for the seeing) 
and carried out aperture photometry for various aperture 
sizes (which allowed us to probe the distribution of the 
flux and the source size in the sky) . We then calculated the 
flux through an artificial 1" slit oriented along the same 
PA as for the observation. The average ratio between the 
total flux and that within the 1" slit gave us an average 
correction factor of ~1.7. 

(2) As a check, we have compared the broad band mag- 
nitudes in the z and J bands (corresponding to the Ha 
rest-wavelength for z~0.8 and z~0.6, respectively) with 
the total observed Ha line flux emitted by the sources. 
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(c) HDF-S: s55. (d) HDF-S: s82. 

Fig. 3. Spectra of the 15 fira sources in the HDF-S observed with ISAAC at the VLT. 



Since we cannot estimate the line-to-continuum ratio from 
our spectra due to the low continuum signal, the ratio of 
the broad-band flux to the Ha flux provides a conservative 
upper limit to the aperture correction. 

The correction factors that we obtained in this way 
are ~ 2 on average for the four sources, consistent with 
the previously estimated average value. In what follows 
we used correction factors of 1.15 for the unresolved and 
1.7 for the resolved sources. 

3.2. Evidence for AGN contributions to the line 
emission 

A primary motivation for our spectroscopic follow up of 
the faint ISO sources was to determine to what extent 
the luminous IR emissions are contributed by AGNs, and 
how much by young stars. To this end, the observed opti- 
cal and IR spectra have been used to identify broad lines 
components and anomalous line flux ratios. 



Our search for broad permitted lines was limited by 
the low spectral resolution of the ISAAC spectra, and 
even more by the low signal-to-noise ratio preventing re- 
liable setting of the underlying continuum. The spectral 
resolving power of R s ~600 corresponds to a rest-frame 
resolution of 12A at z=0.6 and a velocity threshold of 
~ 250 km/s for the ISAAC Low-Resolution gratings. For 
the Medium-Resolution grating the resolving power cor- 
responds to ~ 40 km/s. Consequently, for 16 out of 20 
ISO sources in our sample observed at low-resolution, only 
broad-line components (v> 1000 km/s) can be resolved in 
principle. For S14, S20 and S23 we have good quality op- 
tical NTT/EMMI spectra with resolving power R s = 580 
(~ 500 km/s). 

We report in Table0|the half-power widths of the rele- 
vant permitted lines based on both IR and optical spectra. 
Only for one source (S19) these imply velocity fields signif- 
icantly in excess of 1000 km/s (indicative of the presence 
of an AGN). For four more sources with H a + [Nn] line 
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Fig. [3] continues. 



widths also formally in excess of 1000 km/s, the poor S/N 
and the unresolved [Nil] contribution prevent a reliable 
measure of the velocity field. In all other cases the line 
widths are consistent with those typical of massive star- 
burst galaxies. 

The H a to [Nil] flux ratio may be used in princi- 
ple as an indicator of the ionization field in the source. 
However, because of our limited resolving power, the line 
complex was resolved for only 3 sources of the 1999 run 
(Rigopoulou et al. 2000). Interestingly, one of these (S38 
at z=1.39) reveals an inverted line ratio (H a /[Nii]~0.8), 
together with a moderately broadened H Q (~ 600 km/s, 
see Table El. 

No evidence for broad components or peculiar 
H a /[Nn] flux ratios was found in the optical counterparts 
to S25, S27 and S55 observed in the ISAAC Medium- 
Resolution mode, whose spectra are consistent with those 
of standard spiral and starburst galaxies with no AGN 
signatures (Rigopoulou et al. 2002). 



Altogether, two out of 21 of the faint ISO sources in 
our spectroscopic survey reveal evidence for either type-I 
(source S19) or type-II (source S38) AGN activity. These 
indications from line measurements will be compared in 
Sect 14. 21 with those coming from the study of the optical- 
IR-radio continuum SEDs. 

3.3. Estimates of dust extinction from line ratios 

Corrections of the H a flux for dust extintion are needed 
when using the line fluxes to estimate the rates of star 
formation. These are usually computed from the ob- 
served ratios of the Balmer lines compared with theo- 
retic al models of atomic tran sitions and nebular emis- 
sion. iHummer k. Storevl l)l987|) report the ratios of emis- 
sion lines based on Case B recombination theory for 
T = 10000: 

/(Ha) 



J(H/3) 



2.85. 
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(d) HDF-S: sl9, FORS1 spectrum. 



Fig. 4. Observed optical spectra: panels (a) to (c) show 15 /jm sources observed at NTT/EMMI. Panel (d) represents 
the QSO sl9 from FORS1/VLT: the two bright broad lines are MgnA2798 and Cm]A1908, redshifted at z ~ 1.56 
(courtesy of Rigopoulou et al., 2002, in preparation). 



By comparing this value with t he observed line ra tios and 
using a standard extinction law l|FitzDatrickul999fl . we ob- 
tain the Ay values reported in Table EJ The same Table 
also includes the extinction values based on fits of spec- 
trophotometric models to the optical/near-IR SEDs, as 
explained in Sect. 14. ll 



All Ay values turn out to be greater than 1.5 magni- 
tudes, implying substantial amounts of dust extinction in 
these objects, which are not typical of normal quiescent 
spirals for which Ay — 0.3 — 0.5. If anything, these es- 
timates are likely to be lower limits, since this indicator 
based on Balmer line ratios quickly saturates in a mix of 
dust and radiative sources. 



3.4. Estimates of the Star Formation Rate 

Nebular emission lines, such as H a , are generated in the 
interstellar medium, ionized by the ultraviolet continuum 
of the young stars recently formed in the starburst (those 
with ages < 10 7 yrs). These lines then provide a direct 
indication of the rate of the ongoing stellar formation. 

The conversion factor between ionizing flux and SFR 
may be computed with an evolutionary synthesis model, 
assuming solar abundances, a Salpeter IMF (0.1-100 M ) 
and for continuous bursts with duration in the range be- 
tween fewxlO 7 yrs to fewxlO 8 yrs. We adopt for this the 
Kennicutt (1998) relation: 



SFR 

M Q yr~ 



7.9 -Hr 42 ^ (erg s" 1 ), 



(1) 



which we used to derive the values of ongoing SFR in the 
HDF-S objects with reliable H Q detections. After conver- 
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Table 4. Widths of relevant emission lines from ISAAC 
Low-Res and optical spectroscopy. The source name, ref- 
erence line used for the measure, FWHM in A, FWHM 
in km/s, and the threshold resolving power in km/s are 
reported. 
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T : Ha fluxes from Rieououlou ct al. (2000) 
*: Ha luminosities in units of 10 erg s~ 
x : S73, no aperture correction performed. 



Table 5. Relevant spectroscopic data: aperture-corrected 
Ha luminosities, rates of star formation based on Ha, and 
extinction estimates from the Balmer decrement. In the 
fourth column we report extinction values based on pho- 
tometric fits with spectral synthesis codes (see Sect. 4.1) 
and values of SFR corrected for slit-aperture and redden- 
ing (from Balmer lines when available). SFRs are in Mq 
yr _1 . 



sion from fluxes to luminosities and after applying slit- 
aperture and extinction corrections, SFRs were calculated 
from eq. and reported in Table El 
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* Based on the AGN model described in Appendix A. 



Table 6. Values of the baryonic masses, SFRs and bolo- 
metric IR (8-^1000 (an) luminosities of the 15 /zm sources 
in the HDF-S, as derived from the analysis of their spec- 
tral energy distributions (see Sects. 4.3 and 4.4). 



4. SOURCE PROPERTIES BASED ON THE 
CONTINUUM EMISSION 

The extensive spectral coverage of our source SEDs be- 
tween the UV and the far-infrared (including the radio 
flux for a few) allows adding very significant constraints 
on the physical processes taking place inside them. 

For typical IR galaxies the bulk of stellar formation 
happens inside dust-rich and optically thick molecular 
clouds, absorbing the UV-optical continuum emitted by 
luminous young stars and re-emitting it in the infrared. 
Given the large extinction values, the optical to near- 
IR continuum spectrum is moderately influenced by the 
starburst and includes important contributions by the 
less-extinguished older stellar populations (e.g. Poggianti, 
Bressan, Franceschini 2001). On the contrary, the mid- 
and far-IR spectrum is typically dominated by thermal 
emission by dust mostly illuminated by the young stars 
(the radio flux is also proportional to the number of re- 
cently born massive stars). 

As a consequence, our capability to cover simultane- 
ously the whole SEDs up to the far-IR provides a unique 
opportunity to sample all stellar generations and ages in 
these galaxies, and in particular to compare the mass frac- 
tions in old stars (the inactive, passively evolving compo- 
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nent) to those of newly formed stars (the "active" com- 
ponent), hence to evaluate the "activity" level in these 
galaxies. 

SEDs and (when available) WFPC-2 images of our 
sample sources are reported in Figs. A. 2 and A. 3, together 
with spectral fits based on GRASIL, and with the SEDs 
of template galaxies used for comparison. The dot-dashed 
line in each panels, in particular, corresponds to the SED 
of M51, which we assumed to represent the prototype in- 
active spiral: the comparison of this template, scaled to 
fit the near-IR spectrum, with the ISO data at 6.7 and 
15 /an shows that quite often the observed mid-IR fluxes 
are in excess by a substantial factor. Assuming that this 
excess IR emission is due to the starburst, this factor may 
be taken as a measure of the level of "activity" in the ISO 
galaxies. 

Additional contributions to the source fluxes - fur- 
ther enhancing this "activity" level - may come from 
nuclear non-thermal emission of gravitational origin, an 
AGN component. 



4.1. Estimates of dust extinction based on the 
continuum 

The lack of information on the Balmer line ratios for sev- 
eral of our sample sources has forced us to look for alter- 
native methods to estimate the dust extinction. A widely 
used one is to exploit colour indices. Using the codes of 
spectrophotometric synthesis GRASIL (Silva ct al., 1998) 
and STARS (Sternberg 1998, Thornley et al. 2000) we 
have generated a grid of models for various star forma- 
tion histories and bursts of different duration, and calcu- 
lated the intrinsic (dust-free) V-K colours (typically vary- 
ing in the range V-K = 1.1-1.5). We have then applied 
infrared and optical K-corrections from Poggianti (1997) 
and Coleman (1980), respectively. By comparing the ob- 
served and model predicted V-K colours, we obtained a 
median Ay of 1.8, assuming a screen model for the ex- 
tinction. We stress that the extinction estimates based on 
the V-K colour represent a global (galaxy-wide) extinc- 
tion which may not necessarily represent the extinction 
towards the deeply embedded young stars where most of 
the line emission originates. Note that, although some- 
what more physically motivated, a model with homoge- 
neous mix of dust and stars could not provide a fit to the 
optical spectrum: at increasing the intrinsic dust optical 
depth the extinction saturates to Ay — 0.7 (e.g. Poggianti, 
Brcssan, Franceschini 2001), a value not large enough to 
explain our typical source spectra. 

Alternatively, we can estimate the extinction based 
on the rest-frame Ltjv(2800) continuum. Values of Ltjv 
have been interpolated from B V I colours, using grids of 
synthetic models as described above, with constant star 
formation rate. The median Ltjv (2800) for the present 
ISOHDF-S sample is 1.47xl0 40 erg s" 1 A" 1 . For the same 



IMF as assumed before, we find the following relation be- 
tween the SFR and L uv (2800): 

S F R(M & / yr) = 1.88 x l()- 4a L vv {2800){erg S - 1 A' 1 ). (2) 

The ratio between the H Q and UV-based SFRs turns 
out to be SFR(H Q )/SFR(2800)-4.0. For the correspond- 
ing extinction it follows that A(H a ) = 0.76 Ay and 
A(UV)=1.6 Ay (Pei 1992). We then deduce that, for a 
screen distribution, the extinction is Ay ~1.6. This extinc- 
tion corresponds to a correction factor for the SFR(H Q ) 
of ^3. 



4.2. Disentangling AGN from starburst signatures in 
the UV-optical-IR-radio SEDs 

The mid-IR spectra of galaxies are a complex mixture of 
various components. Typically the most important one is 
the emission by molecular complexes (probably Polycyclic 
Aromatic Hydrocarbons, PAH, see Puget & Leger 1989, 
but their nature is still uncertain) producing bands at 6.2, 
7.7, 8.6 and 11.3 ^m on top of a hot dust continuum. The 
PAH emission is particularly prominent in star-forming 
galaxies. Both the continuum intensity and the PAH's 
equivalent widths are sensitive to the presence of an AGN: 
while the AGN luminous point-like source enhances the 
hot-dust continuum emission, the PAH molecules tend to 
be destroyed by the intense radiative field and the corre- 
sponding emission bands become weaker (Lutz et al. 1998, 
Rigopoulou et al. 1999, Tran et al. 2001). 

Therefore starburst-dominated and AGN-dominated 
galaxies tend to display different mid-IR spectra. The for- 
mer show a bumpy spectrum with prominent emission fea- 
tures at A > 6 /jm and a steeply decresing flux shortwards 
of 6 /im, due to the moderate intensity of the radiative 
field and to the lack of very hot dust. A typical starburst 
spectrum is reported in Figs. A. 2 and A. 3 as the thick 
continuous line, corresponding to the SED of the galaxy 
M82. 

On the contrary, AGN-dominated sources show rather 
flat mid-IR spectra, with almost absent PAH features and 
strong emission by very hot dust detectable down to few 
microns (a comparison of the spectrum of the Seyfert 
galaxy NGC 1068 with those of starbursts is reported in 
Aussel et al. 1999 and Elbaz et al. 2002). 

The shape of the mid-IR spectrum can then in princi- 
ple be used to disentangle between the two power mech- 
anisms. In our case we can exploit the ratio of the LW3 
to LW2 fluxes as a measure of how fast the rest-frame 
SED drops at A < 6 /im for sources at z > 0.4: while this 
flux ratio is expected to be > 4 for starbursts as a conse- 
quence of the LW2 flux missing the redshifted PAH and 
dust emissions, the ratio becomes g^^j — 4 m tnc case 
of an AGN-dominated source. Two such sources in Fig. 



SJLW3) 



~ 1.5 



A. 2 are S19 and S38, having respectively s(LW2) 
and 2.3. For a third source, S82, the flux ratio is also low, 
g^j^j ~ 2.6, and this could also be a type-II quasar. In 
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all other cases, either LW2 has no detection or the flux 
ratio falls in the starburst regime. 

Note that this starburst / AGN discriminant would not 
work for sources at z < 0.4: an example is S73, show- 
ing a low f [^^2) — ^ va hie only because its low redshift 
(z=0.17) prevents the PAH bundle to be redshifted out of 
the LW2 band. 

Radio data could in principle add valuable information 
to the SED analysis. On one side, given the tight radio/far- 
IR correlation for star-forming galaxies (Helou, Soifcr and 
Rowan- Robinson 1985), the radio flux may provide an in- 
dependent estimate of the bolometric luminosity and star 
formation rate. On the other side, a substantial excess of 
radio emission above the value pertaining to the radio/far- 
IR relation may be taken as evidence for emission by an 
AGN. 

Unfortunately, for only a handful of the HDF South 
sources data at 1.4, 2.5, 4.9 and 8.6 GHz are reported 
by Mann et al. (2002) based on ATCA observations. Of 
the four radio-detected sources in our sample, S19 was al- 
ready identified as a bona-fide type-I quasar, while S39 
was suspected to include an AGN contribution based on 
the broadness of the H Q line and morphology (see the 
Appendix). For S39 the 1.4 GHz radio emission is however 
consistent with the radio/far-IR relation for starbursts 
(see Fig. below). 




Fig. 5. Rates of star formation in the 15 /urn sources based 
on the bolometric IR luminosity, as a function of redshift. 
Filled and open circles refer to objects with spectroscopic 
or photometric redshifts. The shaded area represents the 
unobservable region given by our sensitivity limit of 95 
/iJy imposed by the sample selection. 

For the two other radio detections, S23 and S73, the 
radio and mid-IR fluxes are entirely consistent with a pure 



starburst emission (see again Tables [5] and [7| and also 
Mann et al. 2002). 

In summary, our analysis has found clear evidence for 
the presence of AGNs in two of the sample sources (the 
type-I S19, and the type-II S38). For a third source, S39, 
we suspect the presence of an AGN contribution because 
of morphology and a marginally broadened H Q line pro- 
file. For a fourth source, S82, the low LW3/LW2 flux ratio 
could indicate the contribution to the mid-IR spectrum 
of a type-II AGN. For all other HDF South sources any 
AGN contributions should be minor. Although the statis- 
tics is poor, this 10-^20% AGN fraction among the faint 
ISO sources is the same as found by Fadda et al. (2002) 
using the hard X-ray emission as AGN diagnostic. This 
is also in agreement with the results by Alexander et al. 
(2002) exploiting the 1 Msec Chandra X-ray exposure on 
the HDFN, in which 20 of the 41 ISOCAM sources were 
detected, 4 of which were classified as AGNs and 15 as 
emission line galaxies. 

4.3. Estimating the Rate of ongoing Star Formation 

Once the far-infrared activity of the source is reliably at- 
tributed to young stars, a fundamental indicator of the 
source physical status is the rate of ongoing star for- 
mation (SFR). Although UV-optical line and continuum 
emissions are contributed by newly born stars, Sanders & 
Mirabel (1996) have shown that starbursts with bolomet- 
ric luminosities above 1O 1O L0 produce the bulk of their 
energy in the far-IR. For this reason, the bolometric IR 
luminosity of a starburst galaxy is the most direct and 
reliable estimator of star formation. It is then obviously 
important for us to compare our previous estimates of the 
SFR based on the H Q flux with more robust ones based 
on far-IR luminosities. 

To this end, a measure of the source flux around 100 
/jm, where the galaxy IR SEDs are expected to peak, 
would be needed in principle. Unfortunately, this is not 
currently available to us in the HDF South. The imag- 
ing capabilities of ISO are strongly limited at such long 
wavelengths by the poor spatial resolution and high confu- 
sion noise (e.g. Franceschini et al. 2001). There is neither 
much perspective of an improvement until the operation of 
FIRST-Herschel in 2007. An important result by Elbaz et 
al. (2002), however, was to show that the mid-IR flux for a 
large variety of galaxies (from normal galaxies to luminous 
and ultra-luminous dusty starbursts) is extremely well cor- 
related with the bolometric IR emission. They have shown 
that in a large sample of local objects only a very small 
fraction (~ few %) show significant departures from this 
correlation (a remarkable such discrepant case is Arp 220, 
an ultra-luminous IR galaxy [ULIRG] showing a short- 
age of mid-IR flux compared with the far-IR one, due to 
dust self-absorption as evident from the PAH spectrum 
[Rigopoulou et al. 1999, Haas et al. 2001]). 

This almost linear correlation of the bolometric far- 
IR and mid-IR luminosities was proven to hold in local 



14 



Franceschini, Berta, Rigopoulou, et al.: IR Spectroscopy of Faint ISO Galaxies 



1000 



100 = 



x 
3 



0.01 




0.1 - 



1.5 



Fig. 6. Relation of the observed flux in the 15 /im LW3 ISO band versus redshift, as a function of the rate of star 
formation (SFR). This is based on the assumption that the IR flux is dominated by dust-reprocessed emission of 
young stars. The calibration of this relation has been obtained using the M82 spectral template, the detailed response 
function of the LW3 filter (see eqs. 5 and 7 in Franceschini et al. 2001) and the Kennicutt's (1998) relation between 
SFR and L IR for a Salpeter IMF between 0.1 and 100 M Q . 



galaxies by comparing the bolometric fluxes from IRAS 
data with mid-IR fluxes in various channels, including the 
IRAS 12 ^m and the ISO LW3 and LW2 bands: consid- 
ering that for sources at z ~ 1 the rest-frame emission in 
LW2 shifts to the observed LW3 band, this analysis by 
Elbaz et al. proved that the correlation is likely to hold 
for galaxies at least up to z ~ 1 . This is also confirmed by 
the excellent match between the SFR estimates based on 
the mid-IR and radio fluxes for dusty starbursts at z < 1 
(see Sect. O below, and Garrett 2002, Elbaz et al. 2002). 



We applied the Elbaz et al. (2002) prescriptions to es- 
timate the bolometric IR luminosity of our sources. At 
a redshift of z > 0.8 the LW3 band corresponds to the 
LW2 rest-frame, while for z < 0.8 the central wavelength 
of LW3 falls into the IRAS 12 ^m filter. We therefore 
assumed that the luminosities computed from the LW3 
fluxes of the sources in our sample originate from photons 



emitted close to the LW2 and IRAS12 wavebands, and 
then used the Elbaz et al. relations: 

5+2.37 „ i..t \a 7c . n0.998 



L TR = 4.78 H 



x (t>L„[6.75 /H) 
if zaL„ [6.75 \xm] < 5 x 10 9 L 



1.59 



= 4.37" 



-2.35 





x 10~ 6 x (uL v [6.75 urn}) 

9 



1.62 



if jaL„[6.75 /im]>5x 10 y L Q 



vL v \Yb fim] = 0.042 x (uL u [12 urn]) 
Lir ~ 11.1 ■ vL u [15fj,m]. 



(3) 

(4) 
(5) 



From Lm we then computed the SFRs adopting the 
Kennicutt's (1998) calibration: 



= 1.72- 10- 10 Lir [L q ]. 



(6) 



SFR 

M Q yr- 

The results are reported in the fourth column of Table 
|5] Figure[5]shows a plot of these estimated SFRs as a func- 
tion of redshift. The shaded area represents the unobserv- 
able region given by our sensitivity limit of 95 /iJy imposed 
by the sample selection. Error bars have been computed 
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on the basis of the la uncertainties on the LW3 measur- 
ments only, while we did not take into account the scatters 
in the Elbaz's and Kennicutt's relations. The plot shows 
the effect of Malmquist bias induced by the 15 fim flux 
limit of the sample, high redshift objects being detectable 
only if their luminosity, and SFR, are large enough. 

We compare in Figures A. 2 and A. 3 the observed spec- 
tral energy distributions of our sample sources with the 
empirical SED of M82 scaled to fit the measured LW3 
flux (thick continuous line) . The M82 SED has been taken 
partly from Silva et al. (1998), partly from the observed 
ISOC AM-CVF low-resolution spectrum between 5 and 18 
[im bv iFoerster-Schreiber et all l|200lj) . in order to get a 
proper representation of the PAH spectrum which is crit- 
ical for interpreting the LW3 fluxes. In the Elbaz's et 
al. analysis this prototypical starburst falls close to the 
barycenter of the Lfib, to L15 correlation. The SFR values 
appearing in Figs. A. 2 and A. 3 refer to these fits, assuming 
for M82 a SFR ~ 6 M^yr" 1 . 

We report in Fig. H3 our estimated dependence of the 
SFR on the 15 [im flux as measured in the ISOCAM LW3 
filter, as a function of redshift, based on the M82 SED. 
This takes into account in detail the effects of the filter 
transmission function and the complex mid-IR spectrum 
in the K-correction factor. As discussed in Elbaz et al. 
(2002), values of the SFR based on the M82 template tend 
to be lower by ~ 30 — 50% than those based on the more 
sophisticated analysis based on Eqs. [3]to[f)l 

4.4. Estimating the baryonic masses 

A second fundamental parameter characterizing the 
sources of the IR background, that we identified in the 
faint ISO galaxies, is their mass in stars. This integral 
of the past stellar formation activity in the galaxy pro- 
vides a logical complement to the determination of the 
instantaneous rate of star formation. To measure galaxy 
masses in the HDF South, we have perfomed spatially- 
resolved medium-resolution spectroscopy with ISAAC of 
few selected targets. However, the long time integrations 
required have limited such dynamical estimates to 4 galax- 
ies (see Rigopoulou et al. 2002). 

A powerful alternative to the time-expensive spec- 
troscopic investigations makes use of observations of the 
galaxy near-IR SEDs and their moderate dependence on 
the age of the contributing stars (e.g. Franceschini & 
Lonsdale, 2002). This is due to the fact that in a typical 
galaxy the stellar mass is dominated by low-mass stars, 
with evolutionary timescales of the order of the Hubble 
time. As discussed by several authors (Lancon et al. 1999, 
Origlia & Oliva 2000), these moderate-mass stars dur- 
ing the cool giant phase largely contribute to the near-IR 
SEDs of galaxies. However, dusty starbursts show occa- 
sional evidence in the starburst regions for the presence of 
younger red supergiants as shown by the pronounced [CO] 
2.3 /im absorption (e.g. Foerster-Schreiber et al. 2001). 
Obviously, a contribution by such young massive stars, 
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Fig. 7. Dependence of the masses estimated for the HDF- 
S 15 /xm sources on redshift z. Solid circles correspond to 
objects with spectroscopic and open circles with photo- 
metric redshifts. 

even heavily extinguished, would substantially decrease 
the M/L in the near-IR, hence affecting the stellar mass 
estimate. 

We have analysed the UV-optical-NIR data on our 
sample sources by means of a spectral synthesis code 
that we obtained from a modification of that described 
in Poggianti, Bressan & Franceschini (2001) and specif- 
ically devised to model dusty starbursting galaxies. The 
integrated model spectrum has been generated as a com- 
bination of 10 single stellar populations (SSP) of different 
ages: the youngest generations (10 6 , 3 10 6 , 8 10 6 , 10 7 yr) 
and the intermediate (5 10 7 , 10 8 , 3 10 8 , 5 10 8 , 10 9 yr), 
while the oldest populations of stars have been modelled 
as a constant star formation rate for 210 9 <t<1210 9 
yrs. 

The composite spectrum of each SSP made use of 
the Padova isochrones with the Pickles (1998) spectra, 
also complemented with spectra from the Kurucz libraries 
(Bressan et al. 2001). Each stellar generation is born with 
a Salpeter IMF between 0.1 and 100 M© and was as- 
sumed to be extinguished by dust in a uniform screen 
following the standard Galactic extinction law (Rv = 
A V /E[B - V] = 3.1). The extinction value E(B-V) was 
allowed to vary from one population to another, and the 
extinguished SSP spectra were added up to obtain the 
galaxy synthetic spectrum. Twenty parameters in total - 
i.e. the E(B-V) and the stellar mass for each populations 
- are needed to define the synthetic spectrum. Note that 
such large number of independent stellar populations was 
used to get the best possible description of the observed 
spectra and conservative estimates of the uncertainties in 
the model parameters (namely the total stellar mass). 
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This completely free-form, non-parametric model was 
devised to account in the most general way for bursting 
and discontinuous star formation histories characteristic 
of starburst galaxies, as well as for more normal and qui- 
escent formation patterns. It also provides easy implemen- 
tation of different extinction properties for populations at 
different ages. 

We have used this spectral synthesis code to per- 
form an exploratory study of how degenerate are fits to 
the observed UV-optical-near-IR SEDs of IR starbursts 
against variations in the age and extinction of contribut- 
ing stellar populations (Berta et al. 2002, in preparation). 
We have explored the model's parameter space with the 
Adaptive Simulated Annealing method by Ingber (2000), 
including random-number generators, and using \ 2 as a 
goodness-of-fit test. Our simulations have shown that the 
age-extinction degeneracy seriously hamper the mass esti- 
mate, which may be uncertain by factors up to 5 or more 
for some dusty objects. This is partly because the optical- 
near-IR spectral continuum by itself leaves highly undeter- 
mined the contribution of supergiant stars to the near-IR 
flux. 

To better constrain the incidence of young extin- 
guished stellar populations and to evaluate their contri- 
bution to the galaxy's M/L ratio, we have included in our 
spectral fitting procedure also the ISOCAM mid-IR LW3 
flux, which is a good measure of the bolometric emission 
by young stars (see Sect. 14. 3|) . To match the observed mid- 
IR flux, we have computed the far-IR spectrum associ- 
ated with a given solution of our spectral synthesis model 
from the difference between the total non-extinguished 
and the total extinguished UV-optical-NIR flux, and as- 
suming that this absorbed energy is re-radiated in the IR 
with the spectral shape of M82 (see Figs. A. 2 and A. 3). 
At the end, our observable set included the UV, optical, 
near-IR and mid-IR flux data for all sources, as reported in 
Figs. A. 2 and A. 3. Once a best-fitting solution was found, 
the total stellar mass in the galaxy was computed as the 
sum of the contributions by all stellar populations. 

Our estimated values of the baryonic masses M and 
their uncertainties, as obtained with the above procedure, 
are reported in Table El (fifth and sixth columns) and plot- 
ted against redshift in Fig. [7| We see that, although we 
can exploit a good sampling of the galaxy SED, the uncer- 
tainties in the photometric mass estimates are still fairly 
large (typically a factor ^2). Given the very extensive 
exploration of the parameter space that we performed, 
we consider these as conservative estimates of the present 
uncertainties. Tighter constraints could soon be obtained 
from IRAC/SIRTF data between 3 and 10 fim over large 
sky areas, particularly from the SIRTF Legacy Programs 
GOODS and SWIRE (Dickinson 2002, Franceschini & 
Lonsdale, 2002). 

The M values do not show significant correlation with 
z, e.g. compared with the strong observed correlation of 
SFR with z (Fig. [SJ. This reflects our primary selection 
not being on mass but on the SFR value, through the LW3 
flux limit. 



5. DISCUSSION 

5.1. Comparison of independent SFR estimators 

As expected given the nature of these sources selected for 
their strong mid-IR dust emission, the optical estimates 
of the SFRs based on H„ luminosities without extinction 
corrections provide values very significantly smaller than 
those obtained from the mid-IR flux [Fig. |8(a)| . Also the 
H a -based and the IR-based SFR estimates are poorly cor- 
related with each other. 

A large scatter remains even after the He-based esti- 
mates are corrected for line-of-sight extinction, as illus- 
trated in Figure 8(b) (open symbols here refer to Ay esti- 



mates based on fitting of the optical SEDs, filled symbols 
to estimates based on the Balmer line ratios). 
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1 Luminosities at 1.4 GHz in units of 10 22 W Hz -1 . 

2 AGN objects: SFR is not a reliable estimate. 

Table 7. Radio estimate of the SFR of 15 /im sources in 
both HDF's and FF's. In the first column objects named n 
are from the HDF-N sample, those with s in the HDF-S. 
The ISOCAM LW3 sources in HDF-N refer to the cata- 
logue by Aussel et al. (1999). Third column specifies the 
frequency at whi ch the fluxes in the fourth were measur ed. 
Data are from lRichardsl l|2000|) and lMann et all (Hp02). 



An independent test of the SFR can be inferred from 
the radio flux: the radio emission is also unaffected by dust 
in the line-of-sight (though it might be sensitive in prin- 
ciple to free-free absorption at low radio frequencies for 
large column densities of ionized gas, and also sensitive to 
radio-loud AGN components). The relationship between 
the SFR and radio synchrotron emission is established by 
the number o f tvpe-II and t ype-lb supernova esplosions 
per unit time. ICondonl l(l992h finds the relation 



SFR(M > 5M ) L„(1A GHz) 



M Q yri A-IQ 21 W Hz- 1 ' 

For our adopted IMF, this becomes 
SFR(M> O.1M ) L v {lAGHz) 



M Q yri 



1.2 • 10 21 WHz 
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Fig. 8. Comparison between estimates of SFR based on the H Q luminosity and those based on SED fitting of the 
mid-IR flux (as in Table [SJ. In panel (a) the H Q flux is without extinction correction, in panel (b) it is corrected for 
extinction. Lines indicate the relation SFR(Ha) = SFR(IR). In the right panel open squares refer to Ay estimates 
based on fits to the SED, and filled squares to extinction estimated from the Balmer line ratios. Error bars are based 
on flux uncertainties (propagating into the extinction estimate when Ha is corrected). 



To improve the statistics, we consider in this and the 
following sections source samples selected from deep ISO 
observations in both the HDF-S and the HDF-N. The 
ISOCAM LW3 datasets in both areas have been reduced 
in the same way. As for the HD F-N, the ISOCAM data 
reduction bv lAussel et al.1 l|l999f) has detected 77 sources 
at 15 /um, 41 of which have S\^^ m > 100 fiJy and consti- 
tute a complete sample over an area of 25 square arcmin 
(as for HDF-S). HDF-N sources with S 15fim < 100 n.Jy 
do not form a complete sample, but they provide an useful 
extension covering the faint end of the galaxy luminosity 
function. All but 13 of the HDF-N sources have spec- 
troscopic redshifts, while for the remaining we adopt the 
usual photometric estimate. 

Then using published radio fl uxes at 1.4, 2 .5, 4.9 and 
8.5 GHz for HD F-S galaxies bv iMann et all ll2002h an d 
for HDF-N hv lRicha.rds e t"aTl <jl 998|) and iRichard si jiB. 

we computed the source luminosities at 1.4 GHz assum- 
ing for the radio emission a power-law with spectral index 
au, w ith a» — 0-7 (th e mean value for the HDF-N sam- 
ple of I Richards! ^OOfi) ). except for sources S23 and S73, 
for which the spectral index turns out to be 0.5 and 1.0 
respectively. 

The SFR values found in this way, and reported in 
Tableland Fig. El are in good agreement with those based 
on the IR flux. 

In both table and figure we report also the data for 
the sources si 9 and n20 for which we found evidence for 
the presence of an AGN (n20, detected by Chandra in X- 




1 10 100 1000 

SFR IR [M s yr-'] 

Fig. 9. Comparison between the estimates of SFR based 
on the radio flux and those inferred from fits to the mid- 
IR flux. The two methods provide consistent values in all 
cases where the radio is not contaminated by AGN emis- 
sion (as is the case for the two crossed sources). Error bars 
on the IR estimate are based on the flux uncertainties only. 



rays, i s a recognized type-I A GN, see lHornschemeier et alJ 
l|200lj) ; iBrandt et alJ l|200l[) ): a comparison of the radio- 
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based SFR value with those reported in Table HJ1 provides 
an interesting test of an AGN contribution. The good 
match of SFR values based on the radio and the IR con- 
firms the reliability of the latter as a SFR estimator (see 
also Elbaz et al. 2002, Garrett 2002). 

The rates of SF indicated by our analysis for the faint 
ISO sources at z > 0.5 range from few tens to few hun- 
dreds solar masses/yr, i.e. a substantial factor larger than 
found for faint optically-selected galaxies. 

Altogether, our analysis confirms that the mid-IR light 
is a good tracer of the star-formation rate, since it cor- 
relates well with the radio and H Q line fluxes (Figs. |Ht> 
and|nj). On the other end, even after correcting for dust 
extinction, SFR estimates based on the H Q line flux un- 
derestimate the intrinsic SFR of luminous IR galaxies by 
a factor ~ 2, with some large scatter. We have to consider, 
however, that large and uncertain correction factors have 
been applied due to the poor spatial sampling and dust 
extinction effects on the H a measurement, which might 
explain it. Arguments in favour of a large intrinsic scatter 
bewteen optical line and IR bolometric fluxes were dis- 
cussed e.g. in Cram et al. (1998), Poggianti & Wu (2000), 
Rigopoulou et al. (2000), Poggianti et al. (2001), see also 
Goldader et al. (2002). 

This issue will be possibly solved only with the ad- 
vent of the new-generation of optical and near-IR Integral 
Field Spectrographs on large telescopes (e.g. VIMOS and 
SPIFFI/SINFONI on the ESO VLT, see Thatte et al. 
1998), providing data of high spatial resolution and sensi- 
tivity on the optical-UV emissions in such morpologically 
complex systems. 

5.2. Timescales for star formation 

A proper characterization of the evolutionary status of the 
faint ISO source population and their relevance for the 
general process of galaxy formation comes from matching 
the rate of ongoing star-formation SFR with the mass of 
already formed stars. Such comparison is essentially inde- 
pendent of the stellar IMF (the same scaling factor would 
apply to both SFR and M by changing the IMF). 



We report in Figure 10(a) the ratio of the baryonic 
mass to the SFR against redshift for ISO sources in HDF- 
S and HDF-N. For all these we determined masses and 
SFRs as discussed in Sect. 14.41 This activity parameter, 
tsF[yrs] = w^, is a measure of the timescale for the 
formation of stars. The figure shows that, on top of a 
very large scatter, there is an apparent trend for the SF 
timescale tsF to decrease with z (the Spearman rank co- 
efficient is p = —0.44 for the whole sample of 109 sources, 
with a very significant correlation probability [> 99%]). 

On one side, this result indicates that galaxies at z ~ 1 
and larger are more actively forming stars than those in 
the local universe: a lower t$F suggests that the ongoing 
SF is more significantly contributing to the observed stel- 
lar mass, and a larger fraction of stars are being formed 
during the current SF event. Of course, we expect a se- 
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Fig. 10. The timescale of star formation tgp — 
M/SFR [10 9 yrs] of faint ISO sources as a function of 
redshift (panel a), mass in stars (panel b) and SFR (panel 
c). 



lection effect to play, if we consider that at the higher 
redshifts only the IR brightest galaxies can be detected 
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in the flux-limited samples, while the selection is less di- 
rectly influenced by the stellar mass. We expect that our 
main selection bias should operate in preventing detec- 



tion of sources in the upper right corner of Fig. 10(a) less 
active galaxies, those with the highest tsF at the higher 
redshifts, are not detected at 15 /im due to our limited 
sensitivity. 

On the contrary, we do not expect that sources are 
missed in the other corner at lower tsF (lower left side in 
the figure). One effect to consider here is the low redshift 
of the sources and the small sampled volume: although 
detecting a luminous massive galaxy becomes less likely 
here, luminosities and masses suffer a similar bias, so the 
net effect on tsF should be negligible. 



Figures 10(b) 



and 



10(c) plot the star-formation 



timescale against the stellar mass and the SFR. The for- 
mer shows that there is essentially no dependence of t$F 
on M. Some clear segregation is evident in both Figs. 



10(b) and 10(c) between sources brighter and fainter than 
Sis = 100 uJy: the latter are systematically shifted to- 
wards lower values of the mass and SFR, and to higher 
values of t$F (if we exclude a few low-redshift and low- 
mass galaxies). Apparently, the fainter 15 /im sources cor- 
respond to a less "active" class, closer to the quiescent 
spiral galaxy population. 

Altogether, our analysis indicates a trend for a de- 
creased activity of star-formation (per unit stellar mass) 
in galaxies at lower redshifts. 



5.3. Comparison with normal galaxies selected in the K 
band 

We gain further insight into the nature of the 15/im- 
selected population from a match with galaxy samples 
selected in the optical or near-IR. We compare here the 
baryonic masses, while we do not consider the SFRs 
whose estimate based on optical data may be quite un- 
certain. Figures ^] an d [HI compare the distributions of 
stellar masses versus redshift for the IR-selected sources 
with those of K-band selected galaxies with morpho- 
logical classification in the HDF-N and HDF-S (filled 
symbols refer to the ISO galaxies, open symbols to 
the K-selected galaxies). In Fig. ^2 the reference sam- 
ple arc 69 morphologically-classified E/SO galaxies with 
K < 20.15 in the WFPC-2 HDF-N and HDF-S (plus a 
few NIC MOS-HDFS sources), ov er a total area of 11.7 
arcmin l|Rodighiero et all 1200 1|) . 

The comparison in Fig. 1121 is performed with a 
sample of 52 morphologically classified spiral/irregulars 
with K < 20.47 in the 5.7 arcmin 2 WFPC-2 HDF-N 
field l|Rodighiero et all 12000^ . For all datasets, the bary- 
onic masses are consistently estimated from fits to the 
optical/near-IR SEDs, as discussed in Sect. 14.41 (previ- 
ous results by Rodighiero et al. have been corrected for 
consistency with our assumed cosmology). Panels (a) and 
(b) in both figures concern the bright (5i5^ m > 100/1 Jy) 
and faint (Si§^ m < 100/iJy) 15 /xm subsamples respec- 
tively. These figures provide evidence that the IR-selected 
sources include a population of rather massive galaxies 
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when compared with those selected in K, (particularly if 
we consider that, for reasons mentioned in Sect. 14.41 the K- 
flux limit is expected to select the most massive objects). 
This is particularly apparent among the IR brighter sub- 
sample (see panels [a]), whose mass distribution shows a 
tendency to occupy a region of high mass values compared 
with K-band selected E/S0 galaxies. Fainter 15 /im sources 
(panels [b]) appear to be hosted by more moderately mas- 
sive systems. 

Note that, in any case, this comparison of IR-selected 
and K-band selected galaxies has to be taken with some 
care due to the different sizes of the fields covered by the 
ISO surveys (~ 50 arcmin 2 ) and by the near-IR sam- 
ples (5 to 10 arcmin 2 ). These latter arc particularly sub- 
ject to variance effects due to the small sampled volumes. 
The ISO samples themselves show evidence of conspicu- 
ous structure: two overdensities or clusters are evident in 
Figs. EH and E2 at z ~ 0.58 for LW3 sources in HDF-S 
and at z ~ 0.84 in the HDF-N: associated with these we 
can notice galactic mass distributions spread towards high 
values. Mann et al. (2002) report a detailed discussion of 
the cosmic variance and sampling effects in the HDF-S. 

5.4. Nature of the IR-selected galaxy population 

It is evident from the previous discussion that the ISO sur- 
veys above 100 \x Jy select a population of massive galax- 
ies which are associated with active sites of star forma- 
tion. The WFPC-2 images show clear evidence for inter- 



actions and mergers in ^50% of these galaxies (Berta, 
2000; see also Fig. A.2). The redshift survey in the HDF- 
N by Cohen et al. (1999) has also revealed that the LW3 
sources are almost invariably associated with peaks in the 
sample redshift distribution, which may be interpreted as 
either "walls" in the galaxy large scale structure, or galaxy 
groups. In such environment galaxy interactions are max- 
imally efficient. 

How much of these large stellar masses are formed 
during the ongoing event of star formation identified by 
the excess mid-IR emission? Are the currently observed 
episodes of star formation responsible for the bulk of 
the stellar content? Hints can be inferred from the star- 
formation timescales reported in Fig.^J Values of the ac- 
tivity parameter t$F here show a large scatter and rather 
uniform distribution between those characteristic of in- 
active galaxies (tsF ~ I0 10 yrs) and those of very active 
starbursts with t$F < 10 9 yrs. The observed median value 
tsF ~ 1 Gyr, compared with typical starburst timescales 
of ~0.1 Gyr, implies that a single starburst event may ex- 
plain only a fraction of the whole stellar content. A few 
to several such episodes are then required to build up the 
observed large galactic masses, as part of a protracted SF 
history made of a sequence of starbursting episodes on top 
of a lower-level secular SF. 



Another interesting fact illustrated in Fig. 10(a) is that 
the parameter t$F expressing the "activity stage" of the 
ISO sources is rather uniformly distributed between values 
corresponding to maximal activity (tsF £ 10 9 yrs) and 
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those corresponding to almost inactive systems (t$F > 
10 10 yrs). This implies that in between the levels of maxi- 
mal and minimal activity there is a continuity of interme- 
diate stages. 

In conclusion, the mid-IR flux efficiently discriminates 
distant IR galaxies as for not only their "activity level" 
(*sf)' but also for the absolute values of the stellar 
mass and SFR: the brighter (S15 > 100 fJ.Jy) IR-selected 
sources appear to be more extreme in any sense, typically 
more massive, more luminous, and with higher SFR. This 
star formation is likely triggered by strong dynamical in- 
teractions and mergers among the brightest members of 
galaxy groups. High SFRs and already large M values im- 
ply that these are presumably the formation sites of the 
most massive current-day galaxies. 

6. Conclusions 

We have investigated the characters of IR emissions in 
galaxies over a fair cosmic volume between z ~ 0.2 and 1.5. 
In particular, we have systematically exploited the mid-IR 
flux as a best and most reliable tracer of star-formation. 
This study makes use of combined observations with the 
Infrared Space Observatory and ESO VLT of faint mid- 
IR sources selected with the ISOCAM broad-band filter 
between 12 and 18 fim down to Sxs ~ 100 fx Jy and below. 

Previous analyses by Elbaz et al. (2002) and 
Franceschini et al. (2001) have shown that, given their 
areal densities, redshifts, and far-IR spectra, these faint 
sources are very likely responsible for an important frac- 
tion (of the order of 50%) of the energy density contained 
in the Cosmic IR Background between 10 and 1000 /im. 
In consideration of the fact that the CIRB is the major ra- 
diant component in the universe after the CMB and that 
it involves more energy than the UV, optical, and X-ray 
backgrounds put together, we expect that these objects 
have played an important role in the process of galaxy 
formation. 

We have reported on low-resolution spectroscopic ob- 
servations with the near-infrared ISAAC and optical 
FORS spectrographs on VLT for a representative and un- 
biased subsample of 21 objects selected in HDF-S. In ad- 
dition to data on line-emissions and the mid-IR fluxes at 
15 and 7 /im, we have made use of the rich variety of pho- 
tometric observations in the UV-optical and near-IR to 
derive further constraints on the physics of the sources. 
To improve the statistics, we have also used in our anal- 
ysis a sample of mid-IR sources at similar depths in the 
HDF-N observed by Aussel et al. (1999). 

Our main conclusions are hereby summarized. 

— Fairly intense H Q + [NII] emission is detected in virtu- 
ally all the observed sources. The comparison with the 
rip and [Oil], as well as the SED's analysis, indicate 
typically high extinction values Ay ~ 1.5 — 2 to affect 
the line ratios, quite larger than found for local normal 
spirals. While obviously coherent with the IR-selection 
emphasizing excess dust emission, this shows that the 



intrinsic (de-reddened) H Q flux is strong in these ob- 
jects. Enhanced activity is also proven by the fact that 
the mid-IR flux is typically larger by factors > 2 — 3 
than expected for normal galaxies. 

— We have investigated evidence for the presence of 
Active Galactic Nuclei in the core sample of 21 sources, 
as possibly responsible for such enhanced activity, by 
combining all available information: the broadness of 
the Balmer lines, the morphology of the optical coun- 
terparts from the WFPC-2 images, the shape of the 
mid-IR spectrum (from the LW3/LW2 flux ratio), and 
the ratio of the radio to IR flux. We have found un- 
controversial evidence for nuclear activity in only 2 
objects (S19, a type-I quasar at z=1.57, and S38, a 
luminous type-II AGN at z=1.39), while for two other 
objects (S39, a very luminous ULIRG at z=1.27, and 
S82 at z=0.69) we suspect the presence of AGN con- 
tributions. Although the statistics is poor, this result 
showing a low AGN fraction ~ 10 — 20% among the 
faint IR sources is entirely consistent with that found 
by Fadda et al. (2002) based on deep hard X-ray data. 

— Then assuming that for the bulk of these objects a 
starburst is the dominant energy source, we estimated 
the rates of star-formation SFR from all available indi- 
cators (H Q line flux, mid-IR and radio emissions). We 
confirm that the mid-IR light is a good tracer of the 
star-formation rate, since it correlates well with the 
radio and H a line fluxes (Figs. Eb andEJ. We found 
typically high values of SFR^ 10 — 300 M Q /yr for our 
IR-selected galaxies. On the other hand, even after cor- 
recting for dust extinction, SFR estimates based on the 
H Q line show a large scatter and some systematic off- 
set compared with the intrinsic SFR of luminous IR 
galaxies. It still remains to be checked with sensitive 
IR Integral Field spectrographs whether this is due to 
poor spatial sampling by the slit spectrograph, or to 
an intrinsic depletion of the optical emission by large 
dust extinction in the SF regions of the most active 
starbursts. 

— We have exploited the excellent coverage of the SED 
in the optical/near-IR to estimate the second funda- 
mental parameter, the mass in stars M, by making use 
of a newly-devised tool to quantify the M/L of stellar 
populations with different ages and extinction. We find 
that the faint IR sources with fluxes 6*15 > 100 [iJy 
are hosted by massive galaxies (M ~ 10 11 M ), even 
if compared with those selected in the K band at sim- 
ilar redshifts. Spatially-resolved near-IR spectroscopy 
of a few of these galaxies by Rigopoulou et al. (2002) 
supports this result. 

— By matching these large stellar masses with the ob- 
served rates of SF, we have determined timescales for 
SF of tsF ~ 0.2 to 10 Gyrs. When compared with the 
tsF values (~ 10 8 yrs) typically found for starbursts, 
this implies that the ongoing SF can generate only a 
fraction of the stellar content in these galaxies, many 
of such repeated episodes during a protracted SF his- 
tory being required for the whole galactic build-up. A 
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trend towards a reduced level of star-formation activ- 
ity in galaxies at decreasing redshift is also apparent 
in our data. 

In summary, the 15 /urn selection appears to em- 
phasize sites of enhanced star formation inside massive 
galaxies, which are typically the brightest members of 
galaxy groups. These sources probably trace evolution- 
ary phases, involving strong dynamical interactions and 
mergers, bringing to the formation of massive current-day 
galaxies. 
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Appendix A: Notes on individual sources 

Some properties of individual sources are hereby sum- 
marized. Imaging data and SEDs mentioned here refer 
to those of Figs. A. 2. Images are from the HDF South 
WFPC-2 V-606 band observations. 

S14. Outside the area covered by WFPC-2, the [On], 
Hg and [Oin] lines have been detected with EMMI at 
z=0.41. The optical spectrum is quite bluer than that of 
M51, and is fit by the SED of a massive, moderately star- 
forming spiral with SFR=16.5 MQ-yr^ 1 , as also revealed 
by the modest excess of the LW3 15 /im flux above the 
M51 spectrum. 

S16. Also outside the HST fields, it shares very similar 
properties to S14, with an H a detection at z=0.62. 

S19 is an WFP C-2 point-like sourc e, detected also 
at 4.9 and 8.5 GHz iMann et alJ l|2002^ . Its SED (Fig. 
A. 2) is virtually a flat power-law from 0.3 to 20 /im. The 
ISAAC spectrum reveals a strongly broadened (FWHM 
~ 260 A) H Q at z—1.57. The corresponding velocity field 
(Av > 4600 km/s) indicates this source to be a type- 
I quasar. Two broad emission features are also detected 
in the FORS spectrum of Fig. 



4(d) 



corresponding to 
MgnA2798 and Cm]A1908 at redshiftTz ~ 1.56. The ob- 
served SED is reproduced by a combination of a starburst 
template in the optical and a type I AGN model. The lat- 
ter is taken from Granato, Danese & Franceschini (1997) 
and Franceschini et al. (2002), with best-fit parameters 
R = 300 (ratio between inner and outer tori radii) and 
equatorial optical depth t v = 50 mag at 0.3 ^m. 

S20. Very similar to S14 and S16 in all respects, two 
lines detected by EMMI are [On] and R at z=0.39. The 
overall SED is quite well fit by the M51 template, and 
the source looks like a massive but relatively normal and 
quiescent (SFR — 11 M©yr _1 ) spiral galaxy. 

S23 has been observed with NTT/EMMI, FORS2 and 
ISAAC on VLT (ISAAC observations are from Rigopoulou 
et al. 2000). It appears as a triple system in the HST 
image, but unfortunately only one of the three nuclei 
(that in the middle) has been included in the spectro- 
graph slits because of the limited orientation capabilities. 
Four emission lines are detected at redshift z=0.46. The 
central source in the triplet, which we identify as that 
mostly responsible for the IR emission, is extremely red 
(B — K = 5.2). The LW3 flux shows a strong excess and 
indicates SFR = 67 M©2/r -1 . The so urce is also detecte d 
in the radio at 1.4, 2.5 and 4.9 GHz l(Mann et all 12002). 
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Fig. A.l. Spectral Energy Distribution of the object s38. 
The optical-NIR datapoints are fitted by normal spiral 
synthetic spectrum obtained with GRASIL (thin solid 
line, see sect 14. 4[) . while the ISOCAM mid-IR data are 
fitted by a type-II AGN spectral template (dot-dashed 
line) . The latter assumes obscuration by an edge-on torus 
whose ratio between inner and outer radii is 300, and with 
an optical depth at 0.3 /Ltm of r„=30 (see Franceschini et 
al. 2002). 

S25, outside the WFPC-2 frame, is associated with 
a likely interacting pair of galaxies on the EIS images. 
We assumed the optical counterpart to correspond to the 
brighter galaxy. Both optical counterparts have been ob- 
served in the high-res ISAAC mode by Rigopoulou et al. 
(2002), and were found to counter-rotate. H Q is resolved 
from [Nil]. Optical and near-IR spectra consistently in- 
dicate a redshift of 0.58. Our photometrically-estimated 
baryonic mass is ~ 10 Mq, to be compared with the to- 
tal dynamical mass of 4 lO ^M^. 

527 was observed by IRigopoulou et alJ (jlflOO) and 
here with FORS2. H Q and H^g detections indicate z = 0.58. 
Rigopoulou et al. (2002) report results of high resolution 
ISAAC spectroscopy, providing a huge dynamical mass of 
1O 12 M . The optical SED is well fit by the M51 template, 
and yields to a mass of 4.5 ■ lO n M . The WFPC-2 im- 
ages show a prominent bulge and well evident regions of 
star-formation along the spiral arms. The moderately high 
LW3 flux suggests, however, that although this is one of 
the most massive spirals known with evidence of ongoing 
star-formation, overall it looks as a relatively quiescent 
galaxy (see Sect. l4.3|l . 

528 is a substantial, relatively blue, spiral at z = 0.56 
(H Q from IRigopoulou et all [2000) , with a clearly defined 
and asymmetric spiral arm, possibly indicative of an in- 
teraction. 

S38 has been observed bv IRigopoulou et al.1 ((2000) , 
who detected the resolved H a and [Nil] lines. The line 
ratio is inverted ([Nn]A6583 >H Q ), indicating a kind of 
type-II AGN activity. The object is detected at both 
15 /xm and in LW2 at 6.7 /xm: the comparatively large 



value for the latter also clearly indicates an AGN con- 
tribution. If we assume the IR spectrum of a type-II 
quasar (Franceschini et al. 2002) fitting the LW3/LW2 
ratio (see Fig. FOj) . the 8 to 1000 /im luminosity becomes 
Lin = 2.72 • 10 45 erg s -1 . The very red optical source shows 
an almost point-like core with a very faint diffuse exten- 
sion (Fig. A. 2). The moderate broadness of the H Q and 
the red optical colors suggest a type-II AGN nature. 

539 at z=1.27 is the most luminous ULIRG (Ljr = 
4.4 10 12 L©) among the sources in our sample. It is a 
radio source detected at 1.4 GHz. The WFPC-2 mor- 
phology becomes more and mo re point-like at in creas- 
ing wavelengths, from U to K l|Mann et all l2002j) . The 
WFPC-2 image shows a spiral structure with a bright 
nuclear extended emiss ion, a kind of compact bulge. 
IRigopoulou et all l|2000h detect H a at z = 1.27, with 
marginal evidence for broadening (1200 ±450 Km/s). The 
LW3 flux, assumed it is due to a starburst (but an AGN 
contribution is also possible, see Sect. 4.1), leads to a very 
high value of SFR ~ 750 M Q yr-\ 

540 is a strongly disturbed spiral showing three hot 
spots, likely starbursting regions, in the HST image. H a 
is detected in the ISAAC low-resolution spectrum at z = 
1.27. The optical SED is very well fit by the M51 template, 
yielding a large baryonic mass of ~ 2.2-10 11 M Q . The LW3 
flux shows a very large excess, indicative of an actively 
starforming galaxy with SFR ~ 265 M^yr^ 1 . 

S43 is outside the HST field. H a is detected in the 
ISAAC 2000 run at z = 0.95 (the line was not detected by 
Rigopoulou et al. [2000] because of a poor estimate of the 
photometric redshift). The optical spectrum is well fit by 
M51, with a large mid-IR excess. 

S53 is a beautiful double system with clear signs of in- 
teraction. Rigopoulou et al. (2000) published the near-IR 
spectrum of the brighter component. In the optical spectra 
we can deblend two components of Hg, at AA ~ 30 A with 
EMMI and AA ~ 25 A with FORS2. These values yield to 
a difference in radial velocity between the two galaxies of 
about 500 km/s. The brighter galaxy shows a well-defined 
circum-nuclear ring caused by a previous high-velocity 
encounter with the smaller galaxy. This latter shows in 
turn a faint off-center spot on top of a l ow-surface bright- 
ness e mission. The redshift measured bv lRigopoulou et all 
(2000) and confirmed here is z = 0.58. The optical SED is 
rather blue, not well matched by the M51 template (the 
mass estimate is based on a fit with refined spectrophoto- 
metric model). 

S55 is identified with another double interacting op- 
tical source, including a large spiral and a seemingly 
sphero id. H a has been detected by IRigopoulou et alJ 
(2000) in ISAAC low-res and medium-res at z = 0.76, 
and Ha with FORS2 observations. The optical-UV SED 
is resonably well reproduced by M51. The high resolution 
spectrum reveals an "S" -shaped H a emission, and allows a 
dynamical mass estimate consistent with our photometric 
estimate (see Tabic |BJ . 

S60 is the IR counterpart of a high-redshift spiral 
galaxy with a possibly interacting counterpart showing 
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so me evidence of interactio n (a tail). It has been observed 
bv lRigopoulou et afl l)2000|) who detected H Q at z = 1.23. 
The overall spectrum is fairly well fit by the M51 template, 
with only a moderate enhancement of the LW3 flux. The 
galaxy then appears to be moderately active in forming 
stars. However, its photometrically-determined mass ap- 
pears to be enormous (~ 5 ■ 10 u Mq), and similar as such 
to the other system of source S27, apart from the much 
larger redshift. Its luminosity and SFR ~ 265 Mq-l/t^ 1 
are in the ULIRG regime. 

S62 was observed with ISAAC by Rigopoulou et al. 
(2000), who detected H Q at z = 0.73. The source appears 
moderately blue in the optical and quite active in the mid- 
IR. 

572 was observed during the ISAAC 2001 run with 
H a detected at redshift z — 0.55. The galaxy is an edge- 
on spiral in an apparently rather crowded field, but not 
much spatial detail is available in the WFPC-2 image. The 
optical SED is quite red, possibly due to inclination. The 
overall SED is very well fitted by the M82 template. 

573 is a very large nearby spiral galaxy with satellites. 
The H Q line emission is detected in the FORS2 spectrum 
at z = 0.17. The SED has been fitted by us with the 
GRASIL library, the mid-IR spectrum showing excess ac- 
tivity compared with an inactive spiral. The WFPC-2 im- 
age shows a spiral with very well developed arms and a 
bar. Radio e missio n at 1.4 and 2.5 GHz is reported by 
iMann et all l)2002j) . 

S79 H Q line emission at z = 0.74 was detected in 
this source during the ISAAC 2000 run. The optical SED 
matches that of M51, but shows a large excess above it in 
the mid-IR. 

S82 consists of two rather symmetric optical sources, 
one of which almost point-like and the other slightly more 
extended, but with no much evidence of disk-like or oth- 
erwise diffuse emission. Observed with ISAAC in 2000, its 
spectrum yields H Q at z = 0.69. The optical SED is pecu- 
liar, with a sharp change in slope in the I band. The small 
LW3/LW2 ratio may indicate the presence of an AGN in 
one or the other of the optical counterparts. 



This figure "2915fl.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0303223vl 



This figure "fran2915_pag26.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0303223vl 



This figure "fran2915_pag27.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0303223vl 



This figure "fran2915_pag28.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0303223vl 



This figure "fran2915_pag29.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0303223vl 



